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Abstract. Bouncing cosmologies, suggested by String/M-theory, may provide an 
alternative to standard inflation to account for the origin of inhomogeneities in 
our universe. The fundamental question regards the correct way to evolve the 
scalar perturbations through the bounce. In this work, we determine the evolution 
of perturbations and the final spectrum for an arbitrary (spatially flat) bouncing 
cosmology, with the only assumption that the bounce is governed by a single physical 
scale. In particular, we find that the spectrum of the pre-bounce growing mode of 
the Bardeen potential (which is scale-invariant in some limit, and thus compatible 
with observations) survives unaltered in the post-bounce only if the comoving pressure 
perturbation is directly proportional to the Bardeen potential rather than its Laplacian, 
as for any known form of ordinary matter. If some new physics acting at the bounce 
justifies such relation, then bouncing cosmologies are entitled to become a real viable 
alternative for the generation of the observed inhomogeneities. Our treatment also 
includes some class of models with extra-dimensions, whereas we show that bounces 
induced by positive spatial curvature are structurally different from all bounces in 
spatially flat universes, requiring a distinct analysis. 
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1. Introduction 

The occurrence of singularities in classical General Relativity (GR) may be the symptom 
of the existence of a more general physical frame, where GR is embedded as a particular 
low curvature limit. In such a frame, that we will generically call “new physics”, one may 
expect that singularities are cured by the existence of a natural physical cutoff, i.e. a 
scale where the new physics comes into play, replacing GR. Presently, the only consistent 
example of such a general theory is String theory, or its 11-dimensional generalization 
M-theory, where this cutoff is naturally provided by the string length l s . 

Following this reasoning line, in a cosmological context, it is natural to expect that 
the new physics would prevent the Big Bang singularity predicted by GR, allowing 
a continuation of the cosmological solutions before the Big Bang. This idea, firstly 
proposed by Veneziano and Gasperini pQ, has found several alternative implementations 
and recurrently emerges in independent attempts to build consistent cosmological 
models from String-M/theories jEJ GS HI 13- All these models, when formulated in 
the Einstein frame, share the same general behavior: the universe starts from an 
asymptotically flat state and undergoes a contraction that pumps the space-time 
curvature up to the new physics scale. At this point, the new physics masters the 
transition to an expanding phase and the standard Friedman-Robertson-Walker (FRW) 
universe begins. In the present literature, such models are generically called bouncing 
cosmologies. 

The interest in bouncing cosmological models comes from the fact that an 
accelerated contraction seems to solve the flatness and horizon problems as efficiently as 
a standard inflationary phase. Then the question opens whether a bouncing cosmology 
may replace standard inflation once and for all. In order to achieve this fascinating 
objective, a bouncing cosmology should provide its own mechanism of generation of the 
inhomogeneities of the universe, in such a way that all present constraints are satisfied 
at the same accuracy level of standard inflation. The most stringent observational 
constraint comes from CMB anisotropies, which require the initial power spectrum of the 
Bardeen potential (encoding the information on all scalar cosmological perturbations) 
to be (nearly) scale-invariant. Is this obtainable or not in a bouncing cosmology? The 
first studies in the Pre-Big Bang scenario gave a negative answer |^, indicating that the 
inhomogeneities must have an independent origin. For example, they can be generated 
in an axion-dominated era through the curvaton mechanism [2J. 

On the other hand, the authors of the Ekpyrotic/cyclic scenario 5j claimed that 
a scale-invariant spectrum can be generated in a slow contraction era and transmitted 
across the bounce. This has risen a considerable debate [HUH] about the correct way 
to describe the evolution of perturbations through a bounce, but, after several years of 
investigation, it seems quite difficult to draw a general conclusion, valid for all models 
of bouncing cosmologies. 

The evolution of perturbations in a bouncing cosmology can be summarized very 
quickly as follows: in the pre-bounce phase, the Bardeen potential has a growing 
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mode, endowed with a red spectrum (becoming scale-invariant in the limit of very slow 
contraction {3|) and a decaying mode with a blue spectrum (becoming scale-invariant in 
the limit of dust-dominated contraction [ 3 ]). After the bounce, we have two alternatives: 
either the spectrum of the growing mode is inherited by a constant mode in the post¬ 
bounce or the growing mode simply decays in the post-bounce and some other constant 
mode with a more blue spectrum dominates at horizon re-entry. 

If the first possibility occurs, then the Ekpyrotic/cyclic models would become a 
complete alternative to standard inflation. Otherwise, a curvaton mechanism or some 
other one must be implemented to make such models compatible with observations. 
All studies of specific (spatially flat) toy models, where perturbations are explicitly 
calculable, show that the pre-bounce growing mode decays in the post-bounce and 
the Bardeen potential is finally dominated by a blue spectrum, incompatible with 
observations jl01 ITH 12 [1H (the incorrect conclusions reached in Refs. |TH fTH 1 ! are 
revised in Ref. [XU] and reconciled with all other studies; also the model presented in 
Ref. [THj is a particular case of the general model presented in Ref. [THj . with two perfect 
fluids with w a = 1/3 and Wb = 5/3). However, the investigation of specific models does 
not allow to infer any general conclusions, since the possibility that some still unknown 
new physics may intervene at the bounce, producing different effects, remains open. 
Indeed, bounces induced by positive spatial curvature show a different behavior, with a 
generic mixing between growing and decaying modes PUCE]. 

The description of the evolution of scalar perturbations through a bounce has 
acquired an independent interest that goes beyond the Ekpyrotic/cyclic models, thanks 
to the mathematical complexity of the problem, which poses several questions and issues 
of general theoretical interest. 

In this work we present a completely analytical study of scalar perturbations in 
spatially flat bouncing cosmologies, making the only assumption that the bounce is 
governed by a single physical scale. In this way, we are able to establish a univocal 
relation between the final spectrum of perturbations and the form of the total energy- 
momentum tensor. We thus introduce a sort of golden rule for scalar perturbations in 
bouncing cosmologies. 

In Sect. 2 we define the generic bouncing cosmology, stating the hypotheses that 
we need to deduce our general result. In Sect. 3 we elaborate our hypotheses, describing 
the background cosmology. In Sect. 4 we write the linear perturbation equations. In 
Sect. 5 we derive the general form of the effective energy-momentum tensor, including 
the high energy corrections driving the bounce. In Sect. 6 we compute the evolution 
of perturbations in the pre-bounce, in the bounce and in the post-bounce. In Sect. 7 
we illustrate the general results of Sect. 6 by some numerical examples. In Sect. 8 we 
consider possible extensions to models with extra-dimensions and analyze the differences 
between bounces in spatially flat universes and bounces in spatially closed universes. 
Sect. 9 contains the conclusions. 
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Let us start by enumerating our hypotheses: 

(i) We assume that it makes sense to define a four-dimensional metric tensor, both 
in the GR regime and in the new physics regime. Then it is possible to write the 
four-dimensional Einstein equations with an effective energy-momentum tensor on 
the right hand side, encoding all the corrections dictated by the new physics. 

(ii) The universe is homogeneous, isotropic and spatially flat, thus being described by 
the FRW metric. 

(iii) There is a unique physical scale r/B governing the bounce. We will set our units in 
such a way that c— 1, 8ttG = 1, g B = 1. 

(iv) In the GR regime far from the bounce (|//| 1), the energy momentum tensor is 

characterized by a constant equation of state and a constant speed of sound, with 
P/P > -!/3. 

In Sect. 8, we will comment on models with extra-dimensions which may fit into 
hypothesis (i) and discuss the differences with positive curvature bounces. 

Now let us work out these hypotheses. The metric, including scalar perturbations, 

is 

ds 2 = a 2 fg) j(l + 2 cj))dg 2 — 2 B ti dr]dx l — [(1 — + 2 E ti f] dx l dx ] }, (1) 

where we have used the conformal time rj. The scale factor is a(rj), while B, E, if and 
<f> are scalar perturbations. 

The total energy-momentum tensor reads 

™ ( P + $P ~{p + p)v,i \ (2) 

u \ (p + p)V ti -( P + 6p)S ij -i ij )' U 

where 5p is the energy density perturbation, 5p is the pressure perturbation, V is the 
velocity potential and £ is the anisotropic stress. Note that these expressions do not 
require GR, but are just dictated by the symmetry of the problem (hypothesis (ii)). 
Thus, also in the bounce regime, the new physics corrections encoded in T ,L V will 
respect this symmetry and contribute to its components, both at the background and 
at the perturbative level. Since this tensor stays on the right hand side of the Einstein 
equations, it must be compatible with the Bianchi identities, thus it must be covariantly 
conserved: V]fJj = 0. 

The independent background Einstein equations and the energy-momentum 
conservation give the Friedman equations 


m 2 = a 2 p 

( 3 ) 

n 2 + 2 H’ = -a 2 p 

( 4 ) 

p + 3H(p + p) = 0 , 

( 5 ) 
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where the prime denotes derivative w.r.t. 7] and 7i = a!/a is the conformal Hubble rate. 

In a spatially flat universe the Null Energy Condition (NEC), which amounts to 
require p + p > 0, must be necessarily violated during the bounce. Of course, this 
violation cannot be achieved by ordinary matter, but requires the intervention of the new 
physics. One may object that this implies the presence of ghosts at the quantum level. 
While the majority of explicit bouncing cosmologies in the literature indeed contain 
ghosts, there are at least two counterexamples. The Erst is the model by Gasperini, 
Giovannini and Veneziano HH, where the NEC violation is performed by a potential for 
the scalar held involving integration on the spatial volume of the universe. The second 
is the model by Biswas, Mazumdar and Siegel |3j, where the generalized gravitational 
action is an analytic function of the D’Alembert operator acting on the curvature scalar 
without poles. In both cases there are no ghosts in the theory. In any case, the possible 
existence of ghosts for a finite time interval may not be a problem at all na. 

3. Background 

Now let us determine the background evolution far from the bounce, using hypothesis 
(iv). In the pre-bounce we have p = W-p and in the post-bounce p = w + p, with 
constant and w + . The restriction w + > —1/3 excludes the possibility of a post¬ 
bounce inflation, which would spoil the results of the bounce. On the other hand, a 
pre-bounce deflation, excluded by the restriction u>_ > —1/3, would make the setting 
of initial conditions in the asymptotic past problematic. 

We set rj — 0 at the bounce, defined as the time where a'(r /) = 0. The bounce 
duration scale must be of the same order of the curvature scale, according to hypothesis 
(iii). Otherwise, there would be two independent scales. So, the bounce will last from 
7/_ ~ —1 to T] + ~ 1. 

Then, the pre-bounce background functions have the power law behavior 


a{rj) ~ \rj\ q ~ 

( 6 ) 

-tj q- 

V 

( 7 ) 

3 q 2 _ 

Pa - | r? |2+2q_ 

( 8 ) 

q ~ ~ 1 + 3w_ ’ 

( 9 ) 


where we have chosen all constant of integration in such a way that the space-time 
curvature at the onset of the bounce rj ~ — 1 is of order 1. 

In the same way, we can write the post-bounce behavior of the background 
functions, which will be analogous to Eqs. ©-© with all '” in the subscripts replaced 
by 




General solution for scalar perturbations in bouncing cosmologies 

4 . Perturbations equations 


6 


We will write the perturbations equations in terms of the following gauge invariant 
variables 


= if + H{E' - B) (10) 

( = ip + HV (11) 

dp v = 5p- p'V (12) 

dp v = Sp — p'V, (13) 


plus the anisotropic stress £. T is the Bardeen potential, (, dpv and dpv are respectively 
the curvature, energy density and pressure perturbations on comoving hypersurfaces. 
In the following we will drop the subscript “V” from all variables. Then, from now on, 
5p and dp must be intended as comoving energy density and pressure perturbations. 

Using the component (Oi) of the Einstein equations to eliminate 0, the (00), (ij) 
and (ii) components become 


dp = -V 2 T 

q /2 n /! 00/2 o_// 

o .. ri — ri . Zrt / i ___ ___, 

a 2 Hf = -—-C +-—- V I ; + V • 


n 


n 


(14) 

(15) 

(16) 


Some of the coefficients of these equations become singular at the bounce (TL = 0) 
or at the NEC violation time (fH 2 — 7 ~C = 0). In order to check the regularity of our 
variables, we may resort to the regular gauge argument used in Ref. U2t If a gauge 
exists where all components of the metric and energy-momentum tensor stay finite, 
then we can explicitly check the regularity of our gauge invariant variables, expressing 
them in this gauge. As shown in Pi. X I> is regular while ( is not. However, we may 
construct ( = (7 i. 2 — where the extra factor compensates the possible divergence 

of V at the NEC violation, allowed by the fact that it is the combination (p + p)V that 
appears in the energy-momentum tensor and is guaranteed to be regular by the regular 
gauge hypothesis. dp\> has no divergences (since p' — 0 at the NEC violation time, 
thanks to the continuity equation), while for dpv we should use the regularized variable 
dfiv = (7i 2 — 7i')dpv. Dealing with regularized variables is absolutely necessary in all 
numerical applications, because no integrator can step through a divergence. However, 
for our analytical arguments, it is completely equivalent to discuss the equations (ED- 
m or their regularized counterparts. We choose the unregularized variables only 
because they have an immediate physical interpretation, but we would get exactly the 
same results with the regularized variables. 

Going back to our equations, Eq. m represents the analogous of the Poisson 
equation, relating the Lapla.cian of the gravitational potential to the perturbations in 
the energy density (this is not a dynamical equation). Eq. lfTT>l) is a first order equation 
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in £, which depends on the sources 5p and £. The anisotropic stress also enters as a 
source in the first order equation for T (USD- It is easy to verify that the perturbations 
of the continuity equations give no new independent equations. 

All the dynamics of scalar perturbations is contained in these three equations. In 
order to solve them, we need to specify the sources 5p and £ for the energy-momentum 
content of our universe. This is easy for the pre-bounce and post-bounce phases. In 
fact, if they are dominated by a perfect fluid or a scalar field with some potential, the 
sources are 

5p = c 2 s Sp ; £ = 0, (17) 

with s being ” or “+” in the pre-bounce and post-bounce respectively. If the squared 
speed of sound c 2 is equal to w s , we have a perfect fluid, otherwise we have a fluid with 
some intrinsic isocurvature mode fTTT ). In particular, a scalar field with a potential has 
always c 2 = 1, whatever the equation of state induced by the presence of the potential. 
In any case, with these sources, Eqs. (©-(USD can be combined into a second order 
equation for ( or $, suitable for vacuum normalization. Expanding the variables in 
Fourier modes of wavenumber k as usual, in the asymptotic past we have 

c - c ww^^ 1){c - kM) ’ (18) 

where C is a normalization constant of order 1, H^ is the Hankel function of the first 
kind and v = 1 /2 — g_. 

We are interested in modes that are well outside the horizon at the bounce, thus 
having 1. For such modes, in the range —1/k < p < —1, we can use the expansion 
of the Hankel functions for small arguments and Eq. © to determine T. Discarding 
all numerical factors of order one, we have two modes for each variable 

£ ~ k u \r ] \ 1 - 2q - + k~ v (19) 

T ~ k u ~ 2 \rj\~ 1 ~ 2q - + k~ v . (20) 

T is characterized by a growing mode with a red spectrum (becoming scale-invariant 
in the limit g_ —> 0) and a constant mode with a blue spectrum. £ is characterized 
by a mode that is growing for q_ > 1/2 and decaying otherwise. This mode has a 
spectrum that is blue for q_ < 2 and red otherwise, being scale-invariant in the limit of 
dust-dominated pre-bounce j3]. The second mode is a blue constant mode. 

5. General form of the energy-momentum tensor of perturbations 

In order to discuss the evolution of perturbations through the bounce, we need to specify 
the form of the sources. In the ignorance of the new physics governing the bounce, we 
can still select all admissible forms of the pressure perturbation and the anisotropic stress 
and make a complete discussion of what can be expected at the end of the bounce. 

Firstly, let us deal with the pressure perturbation, leaving aside the anisotropic 
stress for the moment. 
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The pressure perturbation can be a generic linear combination of 5p, T, £ and a 
certain number n of additional internal variables 8<pi, with time depending coefficients. 
Of course, the pressure perturbation may also depend on the time and space derivatives 
of these variables. Now, through Eq. (m, any dependence on 5p can be read as an 
additional dependence on T and its derivatives. 

For each internal variable it is necessary to specify an independent linear differential 
equation in order to close the system. We are not interested in any new spectral 
dependence whose origin is different from vacuum normalization of the original scalar 
degree of freedom. Indeed, if such spectral dependences are carried by these internal 
variables, they may play an important role, giving rise to additional isocurvature modes 
in the post-bounce, as in the model discussed in Ref. 0- They depend on the specific 
model and must be carefully considered in any complete treatment. However, our 
question is more general and only involves the fate of the original spectrum of T and 
£. The internal variables may acquire this spectrum only if £, T and their derivatives 
explicitly appear in their equations of motion. But then their feedback on 5p would 
be an effective renormalization of the coefficients of £, T and their derivatives in the 
expression of 5p. 

5p is then just a linear combination of T, £ and their derivatives with time- 
dependent coefficients. However, 5p is scalar under spatial rotations and thus it must 
contain only scalars. The time derivatives of scalars are still scalars under spatial 
rotation, but the spatial derivatives must be combined with themselves or other vectors 
of the theory in order to produce scalars. More specifically, the spatial derivatives can 
only appear as V 2 or A l di, where A 1 is a generic 3-vector of the theory. However, since 
the background is homogeneous and isotropic, all 3-vectors in the background must 
be zero. So, spatial derivatives can only appear in a Laplacian or a function of the 
Laplacian. 

Summing up, the most general expression for Sp is 



where we have supposed that the dependence on the Laplacians is analytic and allows a 
regular power expansion and we have introduced the factor 2/a 2 for later convenience. 
Using the Fourier mode expansion, the Laplacians have been written as —k 2 . With the 
same arguments, a similar expression should also hold for the anisotropic stress 



Of course, any possible crossed dependence between 5p and £ can be easily solved at 
the linear level. 

Then the physics of the model is completely encoded in the time-dependent 
operators F n (rj,d v ), G n (rj,d v ), J n (r/, ££), K n {rp d ri ). Let us make some examples. 

For a single perfect fluid or a scalar field, as mentioned before, the only non¬ 
vanishing operator is G\ = —c 2 . In the bounce induced by two perfect fluids or scalar 
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fields, studied in Ref. na and containing the models by Peter & Pinto-Neto m, 
Finelli ra, and Allen & Wands T2j, two internal variables are present. Taking into 
account their dependence on ( and T, the final expression for the effective energy- 
momentum tensor only contains two non-vanishing coefficients: F 0 , and G i, i.e. the 
pressure perturbation only depends on ( and ArT, while the anisotropic stress is zero. 
The model by Gasperini, Giovannini and Veneziano m is simpler, since it contains no 
internal variables and the only non-vanishing operator is G\ = 4(7 H! — IF 2 )/ip' 2 , where 
ip is the dilaton field in their model. Cartier’s model pm has a very involved effective 
energy-momentum tensor with a non-vanishing anisotropic stress. Using the expressions 
in the appendix of his work, we find that the covariant conservation is fulfilled only in 
the trivial case F = 1/2, of = 0, indicating that some misprint may be present in the 
expressions published therein. So, we renounced to identify our operators F n , G n , J n 
and K n in this otherwise interesting model. 


6. Evolution of perturbations 

In this section we shall study the evolution of scalar perturbations through all 
cosmological eras, making no assumptions on the specific form of the energy-momentum 
tensor, thus retaining all terms in the general expressions m, Ga¬ 
in Ref. [IS] an easy way to find the evolution of perturbations outside the horizon 
was proposed. It amounts to write the first order differential equations in their integral 
form and write the solution as a recursive series. The integral form of Eqs. 03-05) is 


c 


H 

~n2 


c-i(k) 

k 2 


+ 


a 2 (H 2 - H') 


H 2 


c 2 {k) 


a 2 7i 


2{H 2 - n 


-Spdtj 


(dp + J a 2 7i(dp 
r a 2 ?i 


2{n 2 - n 


-k (dp 


(23) 

(24) 


where C\{k ) and c 2 (k) are two integration constants. 

With these equations it is very easy to find the evolution of perturbations in any 
phase. In fact, using the general expressions of 5p and (, we can recast the set of 
equations in the following form 

X(p,k)=X^(p,k) + A(p,k) j B(p,d v ,k)X(p,k)dp, (25) 

where X = (T,(() is the vector of variables, = (ci(k)H/a 2 , c 2 (k)) is the vector 
containing the decoupled solutions and A(p, k ) and B(p, d v , k) are two coupling matrices 
that can be deduced from Eqs. (EH) and (El, once we use the proper form of the sources 
as functions of ( and T. The general solution to Eq. (El is 

oo 

x - = E-y (26) 

2=0 

where AW is determined by A"^' 1 ) as 

X m = A(ri,k) f k)Xi i - 1) (i 1 ,k)dr 1 . 


(27) 
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So, starting from the decoupled solution X(°\ we can construct the full solution of 
the system. Of course, this only works for modes outside the horizon, where we can stop 
to the first few terms of the series because the higher order terms contain extra-factors 
of k 2 -C 1 or k 2 r] 2 -C 1. Let us examine the solution in each phase. 

6.1. Pre-bounce 

In the pre-bounce, the soulution assumes a very simple form. In fact, there we have 
5p = —2c 2 _k 2 ty/a 2 and £ = 0, while the explicit expression for all background functions 
in terms of powers of g is given by Eqs. ©-©• So, it is very easy to perform all 
integrations. Discarding all numerical factors of order one, we start from the decoupled 
solution, which assumes the following form in the power law regime 


^ = k- 2 c 1 [k)\g\- 1 - 2q - (28) 

C (0) =c 2 (fc); (29) 

the first recursion gives 

>P (1) = c 2 (k) (30) 

(31) 

and the second recursion gives 

tf (2) = Cl (k)\g\ l ~ 2q - (32) 

C (2) = k 2 c 2 (k)\g\ 2 . (33) 


We can see that the second recursion gives the same modes already present in the 
decoupled solution, multiplied by k 2 g 2 . Since we are considering the evolution outside 
the horizon, this factor is much smaller than one and we can safely neglect the second 
and all higher order recursions. We are thus left with 

\E r - ~ k- 2 c 1 (k)\ii\- 1 ~ 2q - + c 2 (k) (34) 

C- ~ ci{k)ri l ~ 2q - + c 2 (k). (35) 

Comparing this solution with the expansion of the Hankel functions for small arguments, 
given by Eqs. m and (EDD, we can make the identifications ci(k) = k u , c 2 {k) — k ", 
recovering the two integration constants in terms of the vacuum normalization. 

6.2. Bounce 

Now let us turn to the core of our work. It must be said that the specific shape of 
the solution during the bounce is not our concern, since we just want to know the final 
spectrum of scalar perturbations when they reenter the horizon. However, as we shall 
see in the next subsection, in order to construct the post-bounce solution, we need to 
evaluate the integrals appearing in Eq. (El in a domain covering the whole bounce 
phase, which ranges from g ~ — 1 to g ~ 1 in our normalized units. It turns out that 
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the evaluation of the bounce integrals is very easy, since the k and 77 dependences are 
always factored. In fact, the generic integral has the form 
1 1 

J }{ k , rj)dq = f k ( k ) J fn ( v)dv = fk ( k ). (36) 

-1 -1 

The last equality is crucial. It comes from the fact that the only scales in our 
problem are rj , k and the bounce scale tjb that we have set to 1. Since we have factored 
the k dependence, the only scale left in f v is 77 itself. But integrating over 77 from —1 to 
1 , whatever the dependence of the primitive function on 77, the final result will lose the 
only scale it contained, and would thus be of order 1. It cannot depend on any other 
scale, unless the bounce is actually governed by two or more scales. In that case, we 
cannot solve the integrals in such a trivial way, except in some special limits. 

Since we are going to use several integrals over the bounce phase in the next section, 
we give their explicit derivation in Appendix A. 


6.3. Post-bounce 


In order to build the solution after the bounce, at every recursion we must evaluate 
integrals of the form 
v 


f ( k , Tl ) dri . 


(37) 


— OO 


The integration domain can be split in three parts, covering the pre-bounce, the bounce 
and the post-bounce phase, respectively. We thus have 

T) — 1 1 V 

J f ( k , v ) dr ) — J f ( k , T )) dr )+ J f ( k , rj ) dr ) + J f ( k , rj ) dr ). (38) 

— OO —OO —1 1 

For the pre-bounce integrals, it is sufficient to use the pre-bounce solutions evaluated 
at 77 ~ — 1 . For the bounce integrals, we will exploit the explicit forms calculated in 
Appendix A according to the prescription given in the previous subsection. The post¬ 
bounce integrals can be easily evaluated using the asymptotic background solution in 
the post-bounce phase, where all background functions are replaced by powers of 77. For 
a generic function f(k,r 7 ) = fk(k)\g\ s , we have 
v 


f h(k)\n\‘dn — h(k) 



(39) 


1 

where we have discarded all numerical factors of order 1 , as usual. 

We are now ready to construct the post-bounce solution. Once more, we start from 
the decoupled solution, that in the post-bounce era again assumes the explicit power-law 
form 


= k u - 2 \rj\- 1 - 2q + 

c (0) = k ~ v . 


(40) 

(41) 



-1 
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Then, the first recursion gives 

1 l 

T (1) = \ V \-^ 2g+ k-(\rj\^ + 1) + / f v (v) C (0) + / (42) 

-l -l 

i i 

c (1) = k v {\r]\ 1 - 2q + + 1) + J f v (p)Sp w dr] + k 2 J f v (v)t {0) dv, (43) 

-l -l 

where the bounce integrals can be read from the appendix. 

The second recursion gives 

= \r]\~ 1 ~ 2q+ [k u (\g\ 1+2q+ + [if + 1 ) 
i i \ 

f v (v)dp {0) dg + k 2 j f v (v)^ 0) drij ( \r]\ 1+2q+ + 1) 

i 

fv(v)C { 1 ) dv+ J f v (g)f (1) d7] (44) 

-i 

C ( 2 ) = A .-, + 2 (| ?/ | 1 - 2 g+ + | ?/ | 2 + i ) 

+ k 2 (j f n M c (0) + f dr,! 1 - 2 ’* + 1 ) 

1 1 

+ J fr,{v)dp (1) dr) + k 2 J f v (p)f {1) dig. (45) 

-l -l 

No further recursions are needed, since all new terms are negligible for superhorizon 
perturbations. At the end, according to Eq. (I2HI) we can sum up all contributions to 
build the full post-bounce solution. A careful analysis reveals that almost all terms 
coming from the bounce integrals are of the same order of other terms already present 
or are negligible for perturbations outside the horizon at the bounce. Then, the final 
post-bounce solution takes the amazingly simple form 

T + ~ k u - 2 \rj\- 1 - 2q + + k~ v + k v + G 0 k u ~ 2 (46) 

C+ ~h v ri l - 2q + + k~ u + k u + G 0 k u ~ 2 , (47) 

where the constant Go is set to 1 if the corresponding operator defined in Eq. m is 

non-vanishing; otherwise we have Go = 0. In other words, this term is present only if 

5p is directly proportional to T or its time derivatives, without k 2 factors, during the 
bounce. 


6.J h Discussion 

The whole history of perturbations in a bouncing cosmology is contained in Eqs. (d 
and m for the pre-bounce phase, and (EED, (S3) for the post-bounce phase. 

Let us start our discussion from the pre-bounce. The Bardeen potential in the pre¬ 
bounce is dominated by the growing mode k v ~ 2 \?l\~ 1 ~' 2q -, endowed with a red spectrum 
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becoming scale-invariant in the limit of slow contraction q —» 0. This mode comes 
from the decoupled solution of Eq. (ED, going as TL/a 2 and is destined to decay in the 
post-bounce phase, where it takes the form k v ^ 2 \g\~ l ^ 2q+ in Eq. (ED- There is also 
a subdominant constant mode k~ v , which plays no role in this phase. ( couples to T 
through a Laplacian and thus has a mode going as k v [r]\ 1 ~ 2q+ in Eq. (TTT71) . where the 
potentially scale-invariant spectrum is multiplied by a k 2 factor. This mode dominates 
for a fast contraction q_ > 1/2 and is subdominant w.r.t. the constant mode k~ u , coming 
from the decoupled solution of Eq. ED, in the case of slow contraction < 1/2. 

Now let us examine the post-bounce phase. The role of the bounce is to upset the 
scalar perturbations for a finite time and to determine the initial conditions for the post¬ 
bounce phase. Immediately after the bounce, in Eq. (l-EBl) we have the decoupled mode 
k v ~ 2 | ? /| _1_2y+ , carrying the potentially scale-invariant spectrum, two constants modes k u 
and k~ u that are always present with a blue spectrum, and possibly one more constant 
mode with the k v ~ 2 spectrum. This additional mode is present only if the comoving 
pressure perturbation is proportional to T or its time derivatives without Laplacians. 
Recalling that we are in the limit k ■ Cl, this constant mode will dominate undisturbed 
until horizon re-entry and will determine the spectrum of scalar perturbations seeding 
the inhomogeneities of the universe. 

In the other case (i.e. if 6p is related to T only through Laplacians), the Bardeen 
potential is initially dominated by its decaying mode. But at horizon re-entry, this 
mode has decayed by a factor r]~f ~ 2q+ , where )] re is the time of horizon re-entry. Setting 
7] re ~ k~ l , we have different possibilities. After a slow contraction in the pre-bounce 
phase, the larger constant mode is k~ u . Then, the decaying mode will still dominate 
only if q + < g_, i.e. if the expansion in the post-bounce is slower than the pre-bounce 
contraction. Otherwise, the constant mode dominates. However, in a realistic cosmology 
the relevant modes for CMB enter the horizon in the radiation or dust-dominated era, 
with q + = 1 or 2, respectively. Then, the spectrum would be totally determined by the 
blue constant mode k~ u , that in the very slow contraction limit q_ -Cl gives a spectral 
index n s = 3. The final possibility is that the universe has undergone a fast pre-bounce 
contraction q_ > 1/2, so that the largest constant mode is k v . In this case, the decaying 
mode dominates when the expansion is slow (g + < 1/2), while is negligible otherwise. 

The modes dominating at horizon re-entry in all cases are summarized by the 
following diagram 


Gn — 1 


T ~ k v ~ 2 


Q- < 1/2 

q- > 1/2 


q + < q_ =>• T ~ k u 2 r] 1 2q+ 
q + > q- =>- T ~ k~ v 
q + < 1/2 ^ T ~ k v ~ 2 r]~ 1 ~ 2q+ 
q+> 1/2 ^ T ~ k v 


(48) 


At the same time, we see that there is a direct correspondence between the modes 
present in T and those in (, so that an analogous discussion is valid also for (. It is 
interesting to note that a dependence on T or its time derivatives in the anisotropic stress 
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does not produce any sizeable effects on the scalar perturbations. The anisotropic stress 
enters without Laplacians in the equation for T, where, at most, it may renormalize the 
amplitude of the decaying mode. In the equation for £ it comes through a Laplacian, 
thus being completely ineffective. 

A very important caveat to keep in mind is that the presence of a dependence 
on X I> in 5p is just a necessary condition for the persistence of the k v ~ 2 spectrum in 
a post-bounce constant mode. In fact, bounces with particular symmetries may bring 
to special cancellations in the integrals, erasing the lowest order contributions to the 
constant mode and leaving space to higher order (most blue) contributions. 

Another possibility is that the internal variables of the new physics at the bounce 
may provide new modes characterized by alternative spectra with a different physical 
origin w.r.t. quantum fluctuations in the initial vacuum state. These modes may 
contribute to the constant mode at the end of the bounce, providing another possible 
alternative source for phenomenologically acceptable scalar perturbations. However, in 
such case, the phenomenology would completely rely on the high energy physics at the 
bounce, whereas it would be preferable that the final spectrum is as independent as 
possible of the unknown bounce physics. 

Our conclusions seem quite similar to those reached by Durrer and Vernizzi in 
Ref. i, where the bounce is replaced by a spacelike hypersurface and the evolution of 
perturbations through the bounce is determined by Israel conditions. In fact, they find 
that the persistence of the growing mode of the Bardeen potential in the post-bounce 
era is very general, while its decay occurs only if the bounce hypersurface coincides 
with the uniform energy density hypersurface, which corresponds to force the surface 
tension to be proportional to fc 2 T. Reinterpreting the surface tension as an integral of 
the pressure perturbation during the bounce, this would correspond to tune G 0 = 0 in 
our expression for 5p. Thus, the two approaches seem in complete agreement. 

7. Numeric testing on toy models 

Our golden rule for scalar perturbations through bouncing cosmologies is based on a 
fully analytical derivation. However, it is useful to visualize the main conclusions on 
some toy models of bouncing cosmologies, so that they can be touched “by hand”. In 
order to make things as simple as possible, we have considered a bounce with two perfect 
fluids: the first dominating the pre-bounce and post-bounce eras and the second, being 
characterized by negative energy density, representing the new physics effective tensor 
inducing the bounce. For a bounce with two perfect fluids, the Friedman equations 
are also integrable in terms of hypergeometric functions na, though we will not use 
the analytical solutions. So, in this section, p a is the energy density of the first fluid 
and pb is the energy density of the second fluid. The two fluids have equations of 
state p a = w a p a , pb = WbPb, with w a < Wb, in order to let the secondary fluid become 
negligible far from the bounce. According to Eq. we have g_ = q + = q. A slow 
contraction followed by a slow expansion (q < 1/2) takes place for w a > 1, while a fast 


General solution for scalar perturbations in bouncing cosmologies 


15 


contraction followed by a fast expansion (q > 1 / 2 ) takes place when w a < 1 . As for all 
other models already present in the literature, in this model the post-bounce phase is 
the time-reversal of the pre-bounce phase. We are thus unable to directly investigate a 
fast contraction followed by a slow expansion or a slow contraction followed by a fast 
expansion. Such asymmetric bouncing cosmologies require the implementation of decay 
of the pre-bounce matter and creation of new matter at the bounce, which would make 
the model much more complicated, though not impossible to treat. 

After the background is determined, we have to specify the form of the energy- 
momentum tensor for perturbations. We choose 

bp = -^w a k 2 ty + 3 (w b - w a )(p b + p b )C + G 0 (p b + Pb)^ (49) 

£ =o. a (50) 


In the comoving pressure perturbation 5p we have included three terms. The 
first one is the classical adiabatic perturbation of the first fluid, which is necessary 
to reproduce the ordinary physics in the pre-bounce and post-bounce phases. The other 
two terms are multiplied by the density of the secondary fluid, which becomes negligible 
far from the bounce. The term proportional to ( is the same that one would have by the 
requirement that the perturbations in the secondary fluid vanish in the gauge '0 = 0 . 
This choice is not essential and has been made just for computational reasons. The last 
term is a “new physics” term directly proportional to T. We will switch it on or off by 
setting Go = 1 or Go = 0, respectively. Here we only present a model with vanishing 
anisotropic stress, as this plays no role in the final result, apart from a renormalization 
of some modes. Of course, we have checked that this is actually the case, testing some 
models (not to be presented here) with non-vanishing anisotropic stress at the bounce. 

With this choice for the energy-momentum tensor of perturbations, the equations 
of motion are simply 


C' + (1 + 3w a )nc 


T' + 


m 2 - h' 

H 




: -w a nk 2 ^ + ^ a 2 n(pb + Pb )'b 

w 


(51) 

(52) 


As anticipated before, it is safer to use ( = (7 ~t 2 — H')( in all numerical applications, 
since this variable is guaranteed to be regular throughout the bounce if any regular 
gauge exists. 

Imposing the vacuum normalization as the initial condition, our toy model is well 
defined and we can choose different values of w a , Wb and Go in order to test the 
conclusions reached analytically. The results are summarized in Fig. |T| where we have 
plotted log |'I , | and log |£| as functions of Sign[? 7 ] loga, which is now a common practice 
for the visualization of the behavior of perturbations in bouncing cosmologies mm- 

For Fig. [TJa, we have chosen w a = 1.5 (implying q = 0.36), Wb = 2.5 and G 0 = 0 . 
For Fig. HJ), we have w a = 1/3 (implying q — 1), Wb — 1 and Go = 0. For Fig. Ht, we 
have w a = 1.5, Wb = 2.5 and G 0 = 1. These three cases illustrate almost everything we 
can expect in a bouncing cosmology. The spectral index n s of the perturbations can be 
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Figure 1 . The three possible behaviors of scalar perturbations in bouncing 
cosmologies, (a) 6p ~ k 2 ^, with q < 1/2 (slow contraction and slow re-expansion); (b) 
5p ~ k 2 ^, with q > 1/2 (fast contraction and fast re-expansion); (c) Sp ~ U/, with any 
value of q. 
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Log a 


Figure 2. The same case as Fig. d with a different relation between Sp and T. Here 
we have chosen 5p ~ (pb +pb)e 1 / p6 T. 


evaluated at any time, evolving at least two Fourier modes with different k and fitting 
the power law k na ~ l . The exponent of the ^-dependence is extracted by the logarithmic 
time derivative of the modulus of the perturbation. In this way we can compare the 
analytical predictions for the k and rj dependences with the numerical results in all 
phases. 

For the pre-bounce, we just have to distinguish q < 1/2 (slow contraction) from 
q > 1/2 (fast contraction). In both cases, T is dominated by the growing mode with the 
red spectrum k v ^ 2 \p\~ l ^ 2q -, as can be seen in all three plots in Fig. d ( is dominated 
by the constant mode k~ u in a slow contraction (Fig. d) and by the growing mode 
k v \ri\ l ~‘ 2q - in a fast contraction (Fig. d)- However, when G 0 = 1, we see in Fig. d 
that the pre-bounce constant mode (we are in a slow contraction) is taken over by a 
mode going as k u ~ 2 \q\ a i SO nre time before the bounce. This mode is induced 
by the direct coupling of ( to T that starts to be effective during the approach to the 
bounce. In fact, in the model presented here, even if the term proportional to T in 
Sp is suppressed by (pb + Pb), grows so large that the G 0 - m °de starts to dominate 
already some time before the bounce. Indeed, choosing a function going to zero faster 
than (pb+Pb) far from the bounce, we can delay the growth of f as much as we like. For 
example, choosing (pb + Ph)e 1 / p6 , the growth of ( practically starts right at the bounce, 
but the post-bounce result is essentially unchanged, as shown in Fig. d 

The post-bounce scalar perturbations essentially depend on the presence of a direct 
proportionality between Sp and T. If we set G () — 1, we see from Fig. d that 
immediately after the bounce both Q and T are dominated by a constant mode carrying 
the pre-bounce spectrum of T, i.e. k u ~ 2 . We want to stress that this outcome is 
completely independent of the expansion and contraction rates and of the specific form 
of the background function that appears in the coupling of Sp to T. Choosing (pb + pb), 
(pb +Pb)e 1 / pb or any other function gives exactly the same result, apart from numerical 
factors of order 1. 
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If the Go-mode is completely absent, then the pre-bounce growing mode of the 
Bardeen potential is entirely converted into a decaying mode, leaving no trace in the 
constant modes. Then, at horizon re-entry it may still dominate or it may be negligible 
w.r.t. some constant mode, depending on the contraction and expansion rates. In Fig. 
Hh we have a slow contraction followed by an expansion with the same rate. We see 
that at horizon re-entry the decaying mode of T is of the same order as the constant 
mode of (, in agreement with the diagram (I48j) . where this case is marginal between 
the dominance of k~ u and the dominance of the decaying mode. In Fig. [TJ) we have a 
fast contraction followed by an expansion with the same rate. In agreement with the 
diagram (@BD, the decaying mode is taken over by a constant mode with a spectrum k v . 

8. Remarks on bounces induced by spatial curvature and extra-dimensions 

In hypothesis (ii) we have excluded the presence of spatial curvature. Indeed, this would 
introduce another scale in our problem and in order to be compatible with hypothesis 
(iii), we have two possibilities. 

The first is that the curvature scale is very low compared to the bounce scale. 
Then the universe is practically spatially flat and we can neglect the spatial curvature 
to follow the treatment of this paper. 

The second is that the curvature itself determines the bounce so that the curvature 
scale at the bounce is actually what we have called rjB = 1. Spatial curvature may induce 
a cosmological bounce only if it is positive, because it enters the Friedman equation 
as a perfect fluid with negative energy density and equation of state w = —1/3. The 
remaining source-matter may keep satisfying p+p > 0, since the effective NEC violation 
is performed by the spatial curvature. In order to let the curvature become negligible 
before and after the bounce, we need a > 0, i.e. decelerated contraction and accelerated 
expansion. For example, the pre-bounce and the post-bounce phases may be dominated 
by sources with equations of state w < —1/3. This contradicts our hypothesis (iv) and 
the initial conditions would be problematic without another early transition from an 
accelerated contraction in the far pre-bounce. This possibility was considered by several 
authors d as a toy model for a bouncing cosmology. 

However, even relaxing hypothesis (iv), there is a structural difference between 
a bounce induced by spatial curvature and all other bounces with negligible spatial 
curvature scale. In fact, the Laplace-Beltrami operator in a positively curved space has 
a discrete spectrum k 2 = n{n-\~ 2) (recall that the curvature scale at the bounce has been 
set to 1). Then, all modes are in the situation k > 1, that means that all modes are inside 
the horizon at the bounce. Then our recursive solution of the integral equation would 
require an infinite number of terms and is completely useless. In practice, all modes 
re-enter the horizon during the pre-bounce deflation and oscillate in the approach to the 
bounce. During the inflationary phase after the bounce, they exit the horizon again, 
stopping the oscillations. This behavior has been put in evidence by Deruellc [tS], who 
has described the oscillations analytically by hypergeometric functions. Then, it is not 
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surprising that mode mixing after the oscillations is a generic outcome of such scenarios. 

We have to conclude that bounces with positive curvature have a substantially 
different physics and cannot be used as a step to investigate bounces in spatially flat 
universes. Then it is difficult to believe that they could be useful for any physically 
motivated models. 

ffowever, there is one remarkable feature of spatial curvature that deserves to be 
mentioned. We have already noted that spatial curvature behaves as a perfect fluid with 
equation of state p = —p/3 at the background level. At the linear perturbation level, 
we may try to recast all equations in such a way that they look like Eqs. m-m with 
some effective energy-momentum tensor containing the curvature terms. After some 
manipulation, we find the following form for the perturbations induced by the spatial 
curvature in the comoving gauge 


Sp = 


-T 


dp 
£ = 




1 

a 2 H 2 


(C-n 


CO 


1 

a 2 Ti 


P(T-C) 


(53) 

(54) 

(55) 


Spatial curvature thus represents a very interesting example of some “new physics” 
with a pressure perturbation directly proportional to T and its time derivative without 
the Laplacian. So, even if spatial curvature is not useful for realistic bouncing 
cosmologies by itself, it encourages the conjecture that the new physics driving the 
bounce may contain terms directly proportional to the Bardeen potential in the effective 
energy-momentum tensor of perturbations. 

Perhaps an even more interesting surprise comes out from extra-dimensions. In 
some particular cases, in fact, it is possible to recast D-dimensional held equations in the 
form of 4-dimensional Friedman equations with some effective energy-momentum tensor 
containing the dynamics of the extra-dimensions. One example is the case in which the 
spacetime is factored into a 4-dimensional FRW metric and an internal homogeneous 
manifold, where none of the metric components depends on the coordinates spanning 
the internal space. For example, consider a five-dimensional spacetime with the metric 

ds 2 = a 2 ('g) j(l + 2(j))dif — 2B i drjdx l — [(1 — 2 if)dij + 2 E yi f\ dx l dx 

— & 2 (?/)(l + 2T)dy 2 — 2Wdtdy, (56) 


with the assumption that all scalar perturbations do not depend on the internal 
coordinate y (i.e. we are only including zero modes and neglecting excitations of the 
extra-dimensions). This metric was used by Battefcld, Patil and Brandenberger [20] to 
study bounces with extra-dimensions. It is easy to show that the background equations 
have the form of 4-dimensional Friedman equations with an effective contribution to the 
total energy-momentum tensor coming from the extra-dimension 


Pextra 


-3 Ha~ 2 


V_ 

V 


Pextra & 


-2 





(57) 
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where the (yy) component of the Einstein equations has been used to eliminate T 5 5 . 

In the same way, examining the perturbations equations, we can re-obtain Eqs. 
(Hl-GHD with effective Sp extra , fipextra and fextra containing terms depending on T, ( T 
and their derivatives (W decouples from other perturbations). Thus, T plays the role of 
internal variable of the model. Its dynamical equation is given by the perturbed (yy) 
component, with 5T 5 5 as a source. 

We will not present the complete expressions of the perturbed effective energy- 
momentum tensor here, but the important lesson is that the presence of terms directly 
proportional to T in Sp is generic in this class of models as well. Then, bounces with 
extra-dimensions may be able to transfer the pre-bounce spectrum of T to the post- 
bounce. This is in agreement with the conclusions of Ref. I20j. 

The fact that extra-dimensions perturbations contribute to the pressure 
perturbation was already noted by Giovannini m- Coherently with the behavior 
of scalar perturbations, he also finds that vector perturbations in multi-dimensional 
bouncing universes preserve the red spectrum acquired in the contraction phase, 
in contrast with 4-dimensional bouncing universes, where the red spectrum is only 
transferred to the decaying mode m 

Nevertheless, a complete analysis of scalar perturbations in multi-dimensional 
bouncing universes that also takes into account non-zero modes is still missing. Their 
relevance certainly deserves some investigation, but they cannot fit into the effective 
four-dimensional picture presented in this paper. Further extensions of our methodology 
are certainly possible. 

9. Conclusions 

In this paper, we have presented the general solution for scalar perturbations in bouncing 
cosmologies, i.e. models of the universe where the present expansion is preceded by a 
contraction phase. This solution has been derived solving integral equations where the 
only assumption on the bounce is that it depends on a single physical scale. The “new 
physics” ruling the bounce has been encoded in an effective energy-momentum tensor, 
for which we have made no a priori assumption. 

The final spectrum of the Bardeen potential in the expansion phase depends on 
the presence or absence in the new physics energy-momentum tensor of a direct relation 
between the comoving pressure perturbation and the Bardeen potential, not mediated 
by a Laplacian, as in all known forms of ordinary matter. If such term is present, the 
pre-bounce growing mode of the Bardeen potential persists in the post-bounce constant 
mode and it is possible to have a phenomenologically acceptable scale-invariant spectrum 
in Ekpyrotic-like cosmologies (i.e. models with very slow contraction in the pre-bounce). 
If such term is absent, the pre-bounce growing mode of the Bardeen potential decays 
in the post-bounce and becomes negligible w.r.t. a blue constant mode, unless the 
post-bounce expansion is slower than the pre-bounce contraction. 

The negative conclusion reached by all regular models of 4-dimensional spatially 
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flat bounces studied up to now is due to the absence of such dependence in the energy- 
momentum tensor. Then, in order to support our conclusion, we have presented a 
numerical study of a toy model with dp oc T which transfers the spectrum of the pre¬ 
bounce growing mode of T to the post-bounce growing mode. 

Bounces induced by spatial curvature are structurally different from spatially flat 
bounces, since all modes are inside the horizon for some time interval centered at the 
bounce. Our treatment cannot be applied to this distinct category of bounces, whose 
study is anyway useless for the understanding of bounces in spatially flat cosmologies. 
However, spatial curvature provides a remarkable example of an effective energy- 
momentum tensor where the pressure perturbation is proportional to the Bardeen 
potential without Laplacians. The same key property is also possessed by models 
of bounces with extra-dimensions where the dependence on the coordinates spanning 
the internal space is neglected. These multi-dimensional models can fit into our 
effective four-dimensional scheme. However, before drawing a positive conclusion, a 
more rigorous and complete treatment (taking into account non-zero modes) should be 
envisaged for multi-dimensional bounces. 

It thus seems not impossible to conjecture a form for the new physics driving 
the bounce that realizes the Ekpyrotic dream. The problem now is not in matching 
conditions or any technical issues about the evolution of perturbations, which are now 
definitively settled by this work. The missing ring of the chain is just the existence 
of some explicit model with the required property of the effective energy-momentum 
tensor and that is physically motivated by a consistent high energy theory like string or 
M-theory. When this explicit model is found, then it will be possible to say that a real 
alternative to standard inflation exists. 

On the other hand, a more prudent attitude suggests the possibility that there might 
exist some important physical reasons precluding the possibility that an effective energy- 
momentum tensor of the required form is allowed at the bounce and then the way for 
bouncing cosmologies would be definitively barred. However, it seems very hard to us to 
find any physical reason to forbid an effective pressure perturbation directly proportional 
to the Bardeen potential. Moreover, the explicit examples of spatial curvature and 
extra-dimensions seem to indicate that this possibility is actually open for high energy 
theories. 
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The calculation of all bounce integrals according to the prescription of Sect. I(i. 2 l is 
straightforward, though a bit boring. Though we need to calculate a relatively large 
number of terms for the second recursion, luckily most of them are similar to other 
terms coming from the pre-bounce or post-bounce phases, so that they are irrelevant 
for the final post-bounce spectrum. 

The starting point is the decoupled solution 

T ( °) = k v ~ (A.l) 

a 2 

C (0) = k~\ (A. 2 ) 

Note that we do not specify the shape of a and TL during the bounce phase. Then, 
indicating by a generic function of 77 , the integrals involved in the first recursion 

are 

1 

//,MC <0) *) = k- (A.3) 

-1 

\ OO 

/ Uv^dv = £ (Jnk~ v+2n + K n k"~ 2+2n ) (A.4) 

_1 71=0 

OO 

/ Uv)dp {0) d V = £ (F n k~" +2n + G n k "- 2+2n ) , (A.5) 

_ 1 n=0 


where the first and the second integrals enter while <£) needs the second (multiplied 
by k 2 ) and the third. We have used the expressions (EHT) and (j221l for bp^ and By 
the symbols F n , G n , J n , K n we generically indicate the results of all integrals involving 
the corresponding operators defined in Eqs. m and (J22D- Of course, these integrals, if 
present, are of order 1, according to the argument in Eq. (ESI). We will keep the symbols 
F n , G n , J n , K n just to remember that these terms are present only if the corresponding 
operators are not zero. At the end it will be sufficient to put these symbols to 1 or 0 to 
consider models where these operators are present or absent, respectively. 

In the second recursion we need the integrals 

\ OO 

/ UvK^dv = £ (F n k-" +2n + G n k"~ 2+2n + J n k ~ v+2n+2 + K n k" +2n ) (A. 6 ) 

_1 71=0 

\ OO OO 

J um {1) dri = E Knk~ u+2n + E (■ JmF n k~ u+2m+2n + J m G n k"~ 2+2m+2n 

1 n= 0 n,m= 0 


I 7 t F,—v+2+2m+2n . t ts FV+2m+2n 
i o i J 777,-TV^a/ 

+K m J n k~ v+2m+2n + K m K n k u ~ 2+2m+2n ) 


(A.7) 


« OO OO 

J fMSp^drt = £ G n j~ u+2n + £ (F m F n k-‘' +2m+2n + F m G n k"~ 2+2m+2n 

n =0 n,m =0 
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+F m J n k- u+2+2m+2n + F m K n k u+2m+2n 

+G m J n k~ v+2m+2n + GmK n V~ 2+2m+2n ) . (A. 

These integrals are those used in Sect. 16.31 to derive the post-bounce solution for 
scalar perturbations. 
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